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I.  INTRODUCTION 


This  Is  the  second  quarterly  report  for  Contract  No.  N000 1  ^4 -8 1  - C -06 15, 
"Application  of  Josephson  Junction  SQUIDs  ^nd  Arrays,"  covering  the  period 
from  2  January  1982  through  2  April  1982.  The  objectives  of  this  program 
are  to  define  and  demonstrate  applications  of  SQUIDs  and  SQUID  arrays  via 
analysis  and  experiment.  The  goals  of  this  contract  are  to  investigate 
and  demonstrate  the  properties  of  a  SQUID  vol tage-controlled-osc i llator 
(VCO) .  These  goals  are  being  met  by  1)  analysis  and  simulation  of  a 
voltage-clamped  resistive  SQUID  and  a  voltage-clamped  dc  SQUID,  2)  design, 
fabrication  and  measurement  of  appropriate  microwave  matching  transformer, 
and  3)  design,  fabrication  and  measurement  of  a  SQUID  VCO.  We  report  on 
progress  in  these  three  areas. 


2. 


ANALYSIS  OF  VOLTAGE-CLAMPED  SQUIDS 


We  have  Investigated  the  voltage-clamped  SQUID  in  two  forms:  the 
resistive  SQUID  of  Figure  1  and  the  dual  resistive  SQUID  with  common 
voltage-biasing  resistor  which  we  call  the  voltage-clamped  dc  SQUID 
(Figure  2).  To  the  best  of  our  knowledge,  the  resistive  SQUID  incor¬ 
porating  the  full  RSJ  model  with  appropriate  capacitance  has  not  been 
previously  investigated  and  this  was  undertaken  to  predict  the  operating 
characteristics  before  final  design  and  fabrication. 

The  circuit  of  Figure  I  is  driven  by  a  current  I  which  produces  a 
dc  voltage  across  r  when  |l|>ic.  The  Induced  oscillation  produces  a 
time-dependent  voltage  across  the  L,C,i  ,R  combination.  We  envision  R 
to  be  the  load  to  which  power  is  delivered.  Hence,  a  determination  of 
this  voltage  will  determine  the  power  delivered.  This  circuit  obeys  the 
relation 

0+(nR+nr)0  +(l+nRnr+ecosQ)0+nf.Ssin0=31+rirP,  (1) 

where  0  is  the  junction  phase 
nR  =  A7C  /R 

nr  =  r/  /L/C 

B  *=  2irLic/<J>o 

gj  =  2ttL  I  /4>o 

and  time  is  measured  in  t/uB  where  ue  =(LC)  .  The  voltage  across  the 
junction  and  load  resistor  R  is  related  to  0  by  the  Josephson  relation 

V=wo<{io0/2ir  .  (2) 

The  last  term  on  the  right  hand  side  of  Eq.  (1),  nrf3  j  =2nr  I  /<f>0w0 ,  is 
approximately  the  expected  Josephson  frequency  <o  ,Ao0  for  l>>i  .  This  is 
actually  equivalent  to  a  voltage  source  (rl)  in  series  with  r;  the  actual 
dc  voltage  across  the  junction  is  somewhat  lower  than  rl  and  must  be 
determined  by  either  direct  measurement  or  calculation  of  <Q>. 


The  simulations  conducted  here  were  carried  out  for  selected  values 
of  o  .  nD,  3  and  6,  with  B.=0.  A  number  of  interesting  and  unexpected 
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results  were  obtained  in  numerical  computation  of  0(t)  and  0(t)  and  are 

reported  here.  The  results  are  essentially  independent  of  n  for  all 

-3  r 

reasonable  values  of  interest  and  so  n  =10  is  used  in  all  calculations 
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reported  here.  This  means  that  I  =  10  i  and,  hence,  6,  =  10  at  u  .=u>©  • 

^  |  C  I  J 

We  also  note  that  hR=Q  ,  the  loaded  circuit  quality  factor. 

The  transient  response  is  dominated  by  hj0,  Q  and  Figure  3  shows 

the  time  evolution  of  0  and  0  for  0=1,  3  =1  and  expected  Josephson 

frequencies  of  uio/20  and  wa .  In  each  case  the  turn-on  transient  is  short 

(approximately  Q./w0)  and  associated  with  shock  excitation  of  the  resonant 

LC  circuit,  and  the  fundamental  frequency  of  the  periodic  response  is  =  w,. 

-1  J 

The  oscillation  shape  appears  to  be  dominated  by  w0  .  We  point  out  that 
the  periodicity  in  a  smoothly  stepped  0  is  2it.  If  the  0  is  raised  to  5, 
the  transient  response  is  longer  and  the  resonance  frequency  becomes  further 
involved  with  the  oscillation.  Figure  4  shows  the  response  at  0=5,  3  =1 
and  expected  ojj /w0  -  1.2  and  1.5  For  Wj/uo  =  1.2  (3 ( =1200),  the  response 
is  anharmonic  although  clearly  there  is  a  strong  component  near  Uj/2we 
and  Wj/wo.  At  u>j/w<,~1.5  ( 3 | = 1 500) ,  the  response  is  periodic  again  with 
a  strong  subharmonic  component  evident  in  the  time  response.  At  (Dj/wo^ 
($1=2000),  the  response  appears  to  be  totally  at  w0=Wj/2,  while  for  Q=1 
there  is  no  evident  subharmonic  response  as  shown  in  Figure  5.  For  those 
situations  where  the  response  is  aperiodic  or  has  strong  components  at 
other  than  ojj,  the  time-dependence  was  computed  over  many  hundreds  of 
cycles  to  eliminate  any  transients  in  assessing  the  oscillator  performance. 
Generally,  for  either  large  3  or  large  Q,  the  response  becomes  very  complex 
and  even  aperiodic.  The  values  required  are  not  I00's,  but  of  the  order  of 
5.  Even  for  smaller  values,  peculiarities  occur  somewhere  in  the  spectrum 
if  3*Q>2. 

In  order  to  classify  and  quantify  these  measurements  we  have  computed 

the  power  spectrum  for  selected  values  of  Q  and  S  as  a  function  of  3(. 

These  were  computed  over  100  periods  of  the  fundamental  frequency  of  the 

oscillator  after  a  long  time  lapse  from  turn  on,  approximately  10^  cycles. 

• 

The  Fourier  amplitudes,  including  the  zero  frequency  or  dc  value  of  0, 
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were  calculated.  The  zero  frequency  amplitude  was  then  compared  to  the 

fundamental  frequency  which  was  first  calculated  in  establishing  the  time 

interval  for  the  numerical  integration.  In  those  cases  for  which  the 

fundamental  frequency  is  the  expected  Josephson  frequency,  <0>=w  .  Other- 
•  * 
wise  <0>=!Dj/n,  except  for  anharmonic  cases  for  which  no  fundamental 

frequency  was  found.  Since  the  Fourier  analysis  was  performed  by  a  Fourier 

series  rather  than  a  Fourier  integral,  no  spectra  were  determined  for 

anharmonic  cases. 

The  power  delivered  to  the  load  resistance  R  at  any  frequency  is 

P,  “  V?/R  (3) 

where  Vj  is  the  Fourier  amplitude  of  the  voltage  at  that  frequency.  From 
Eq.  (2), 


P  .  (5) 

1  (2n)2QL 

•  2 

The  values  directly  computed  are  (0./Q)  which  must  be  multiplied  by 

2  2  * 

[w0<|>o  /  (2tt)  L]  to  obtain  power  in  physical  units. 

Figure  6  shows  the  dependence  of  the  power  delivered  at  Uj  for 

<u=Wo  as  a  function  of  Q  for  8=1 .  Notice  that  the  power  has  a  maximum 

J  -1 
near  Q=3.  The  abrupt  decrease  above  Q=3  is  presumably  the  Q  dependence 

of  Eq.  (5).  The  decrease  at  low  Q  reflects  the  loading  of  the  SQUID  and 

a  strong  decrease  in  0.  The  maximum  value  of  the  power  is  approximately 

2  2  10  - 1 2 
1.124  a><,<f>0/(2ir)  L.  For  design  values  u>©=2  x9x!0  and  1=5x10  H,  this 

corresponds  to  12.9nW.  However,  as  we  have  seen,  values  of  Q>2  leads  to 

some  undesirable  spectral  response.  Nevertheless,  even  for  Q*l,  the 

power  is  not  even  diminished  by  a  factor  of  2  from  the  peak  value. 

Figure  7  shows  a  compilation  of  the  computed  powers  of  the  Josephson 
oscillations  as  a  function  of  Q  for  8=1.  In  the  lower  Q  regime  up  to  2, 
the  spectra  are  well  behaved  except  for  a  subharmonic  response  at 
Wj=l  •  2too  |  Q»2.  For  Q>2,  there  is  a  very  strong  subharmonic  response  at 


Wj*2u)0,  and  anharmonic  response  near  Wj=l .2w0 ,  Q=5.  At  2<Q<5,  the  region 
near  o>j®1.2u)0  exhibits  a  subharmonic  behavior  which  is  not  plotted.  The 
long  dashed  line  with  negative  slope  illustrates  the  approximate  rolloff 
of  the  spectral  harmonic  amplitudes,  roughly  one  decade  in  power  per 
octave.  The  harmonic  response  is  not  plotted  for  clarity. 

The  strong  subharmonic  response  when  <^=20),  which  occurs  at  such 
modest  Q's  as  2.5  suggests  degenerate  parametric  oscillations  such  as  those 
predicted  in  the  SQUID  parametric  amplifier  [Silver,  et  al.,  "Parametric 
Properties  of  SQUID  Lattice  Arrays,"  IEEE  Trans.  Magn.,  MAG-17,  pp.  *t  1 2- 
**15  (1981)].  This  is  also  suggestive  of  the  period-doubling  effects 
attributed  to  strange  attractors  which  are  reported  to  visit  nonlinear 
systems  such  as  Josephson  junctions.  However,  we  also  observe  other  sub¬ 
harmonic  responses  for  Q=5,  8=1  which  do  not  fall  in  the  period-doubling 
class.  An  example  of  this  is  shown  in  Figure  8  for  Q=5,  3=1 ,  and 

8j=l 100. 

These  non- period-doubling  subharmonics  and  the  nonharmonic  response 
which  occur  in  the  general  vicinity  of  Wj=1.2u<>  for  Q=5,  6=1  are  without 
explanation  at  this  time.  A  possible  clue  to  this  behavior  may  be  in 

noting  that  the  resonance  frequency  of  the  circuit  which  we  have  called 
-h 

wo=(LC)  is  in  fact  a  time-dependent  quantity  through  the  phase-dependent 
Josephson  inductance.  The  effective  circuit  inductance  is 

=  L  *  (l+8cos0)  .  (6) 


and  the  resonance  frequency  varies  as 


<*>r  =  (Do  (l+8cosQ)^  . 


Over  each  period  as  0=O-»-2ir,  varies  from  for  6=1.  Since  the 

values  near  u>r=0  occur  when  Q=it,  the  SQUID  spends  a  relatively  short 
time  at  these  very  low  resonance  frequencies  and  a  correspondingly  long 
time  near  «j  Perhaps  some  interactions  between  the  Josephson 

oscillations  and  the  time-dependent  resonant  circuit  induces  the  unexpected 
subharmonics  and  anharmonic  behavior.  We  hope  to  explore  ithis  effect 
further  during  the  next  year. 
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The  effect  of  Increasing  0  at  constant  Q=1  is  illustrated  in  the 
power  spectra  of  Figure  9.  The  bold  printed  points  are  the  Josephson 
frequencies  and  the  lighter  points  connected  by  straight  lines  are  the 
predicted  harmonics.  In  order  to  preserve  some  clarity,  not  all 
harmonic  spectra  are  shown.  Below  3=3. 1A,  all  computed  points  are  well 
behaved  with  some  expected  harmonics  but  no  subharmonic  response.  For 
3=3.1^,  we  observe  subharmonics  as  we  did  for  large  Q  values  at  8=1 . 

At  even  larger  8-values,  the  time  dependent  response  becomes  even  more 
complex  and  we  have  not  computed  spectral  powers.  Such  8-values  would 
be  considered  outside  the  normal  range  for  a  SQUID.  We  report  these 
effects  and  point  out  that  increasing  8  moves  toward  the  non-SQUID 
Josephson  junction  for  which  this  analysis  indicates  very  poor  quality 
VCO. 
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3.  VOLTAGE-CLAMPED  DC  SQUID 


We  have  carried  out  some  preliminary  simulations  of  the  voltage- 
clamped  dc  SQUID  as  shown  in  Figure  2.  The  two  coupled  equations  which 
describe  this  device  are 

0++  (hR+2nr)Q  +  +  [l+2nRnr  +  j  cos  (0++0_)  +  y  cos  (0+-0_)]0  +  + 

[|-  cos  (0++0_)  -  |  cos  (0+-0_)3  0_  +  nr6sin  (0++0_)  +nr&sin  (0+-0_)  =  nr6«  (8) 


0_  +  nR0_  +  y  sin  (0++0_)  -  |sin  (0+-O_)  +  0_  =  0j 


V°L 


V°L 


where  0O  =  2ttL 1 0/4»o  ,  8j  =  2irL  1  j /<}> 0  ,  0+  =  — ^ —  and  0_  =  — ^ —  • 

0R  and  0^  refer  to  the  phases  of  the  right  and  left  hand  junctions, 
respectively,  and  we  have  defined  0  as  2rr L i c/ ({> o  ,  where  L  is  the  inductance 
of  one-half  of  the  dc  SQUID.  In- a  limited  number  of  simulations  with  6=1 
and  6|=0,  the  dc  SQUID  behaves  identically  to  the  resistive  rf  SQUID  with 
0=0=0.  The  two  junctions  switch  in-phase  with  one  another  and  the 
resulting  currents  in  L  add  in  the  small  resistance  r.  For  0j=ir/2  and 
0=1,  the  computation  shows  that  0_  responds  at  the  expected  Josephson 
frequency  while  0+  responds  at  the  second  harmonic.  This  frequency 
doubling  is  a  result  of  the  alternate  switching  of  the  two  junctions  in 
each  Josephson  period.  Figure  10  shows  the  time  dependence  of  the 
various  phases  for  nr=10  nR=QR'=l,  6=1,  0j=it/2  and  0o=IO^. 
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k.  VCO  DESIGN 

The  design  of  the  VCO  circuits  have  been  completed.  The  resistive 
SQUID  is  placed  at  the  low  impedance  end  of  the  500  to  10  transformer 
previously  designed  and  reported.  Connection  to  the  transformer  is  with 
10  microstrip.  Two  devices  will  be  fabricated  on  each  1  cm  x  2  cm 
silicon  chip,  with  k  chips  on  each  2"  wafer. 

Because  of  the  predicted  response  discussed  in  Section  2,  four 

different  designs  are  employed  to  test  the  effects  of  Q  and  L  on  the 

power  output  and  frequency  response.  The  basic  design  is  to  choose  L 

2 

in  order  to  set  the  flux  quantum  energy,  4>0  /2L,  at  a  sufficiently  large 

2  - 1 

value.  Then  set  3=2ttL i c/*4>o=l  ,  which  fixes  ic>  Setting  oj0=(LC)  to  the 
operating  frequency  tunes  the  SQUID  and  selects  C.  If  the  Josephson 
plasma  resonance  is  set  equal  to  to*,  then  the  Josephson  current  density 
is  chosen.  For  the  VCO,  we  assume  that  the  load  resistance  dominates 
over  any  internal  damping.  With  the  previously  designed  transformer, 
this  value  should  be  In. 

The  initial  design  for  the  VCO  was  L=5pH,  £M  ,  <o0=2ttx9x1o'®Hz ,  which 
for  a  in  load  produces  a  loaded  Q=3.3-  Because  of  the  predicted  effects 
for  Q>2,  we  have  designed  four  VCO's  as  follows: 

VCO  1  L=5pH  ,  C=65pF  ,  Q=3.3  ,  ic=67pA 

VCO  2  L=10pH  ,  2C  ,  Q=1.7  ,  2i£ 

VCO  3  L=5pH ,  C  ,  Q=1 .7  ,  !c 

load  with  additional  in  resistance 

VCO  k  L=2 . 5pH ,  2C  ,  Q=1.7  ,  2ic 

load  with  additional  0.25n  resistance 
Although  a  substantial  amount  of  microwave  power  will  be  lost  in  the 
additional  damping  resistances,  the  last  two  VCO's  will  be  devices  damped 
to  a  low  Q  over  essentially  all  frequencies,  whereas  the  damping  provided 
by  the  transformer  in  designs  1  and  2  is  only  well  defined  over  a  few  GHz 
near  resonance.  A  comparison  of  2  and  3  will  be  particularly  instructive. 
Figure  11  shows  a  composite  printout  of  the  seven  masks  designed  for  the 
VCO  wafer. 
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5.  IMPEDANCE  TRANSFORMER 

Work  has  continued  on  the  impedance  and  dimensional  transformer  to 
couple  the  SQUID  to  a  50ft  line.  In  the  last  report,  we  reported  on  the 
electrical  and  geometrical  design  of  the  transformer,  including  the 
computer  generated  prediction  of  the  VSWR  and  a  composite  drawing  of  three 
different  test  transformers.  Ouring  this  reporting  period,  the  photo¬ 
lithographic  masks  were  produced  and  a  number  of  test  transformers 
fabricated  and  measured. 

The  transformers  were  fabricated  in  a  1  cm  x  2  cm  format  on  1" 
diameter  silicon  wafers  with  an  insulating  SiO^  surface.  Three  test 
chips  are  produced  on  each  wafer.  Following  fabrication,  the  wafer  is 
cleaved  with  a  modified  Tempress  scriber  and  the  chips  mounted  in  a 
holder  suitable  for  connection  to  an  OSM  coaxial  bulkhead  launcher. 
Measurements  are  made  at  T=kK  with  an  HP-8klO  Network  Analyzer.  The 
data  are  in  the  form  of  the  scattering  matrix  coefficients  S ^  and 
Only  the  magnitude  of  S  is  reported  since  the  phase  is  difficult  to 
calibrate  with  long  coaxial  lines  into  the  helium  dewar. 

Figure  5  of  the  last  quarterly  report  (15  January  1982)  shows  three 
different  test  transformer  chips.  These  are  mounted  with  a  coaxial 
connector  at  each  end  as  shown  in  Figure  12.  As  anticipated,  connecting 
the  coplanar  input  to  the  bulkhead  connector  was  a  significant  problem. 

An  acceptable  solution  is  to  use  In  solder  to  connect  gold  ribbon  from 
the  connector  centerpin  and  bulkhead  to  the  coplanar  lines.  To  facilitate 
soldering  to  the  Nb  film,  a  gold  film  is  evaporated  over  the  Nb.  This  is 
done  before  the  Nb  is  patterned.  Figure  13  shows  the  response  of  the 
simple  tapered  coplanar  lines  at  k  kelvin;  Figure  lk  the  response  of  the 
double-ended  transformer,  in  which  the  50ft  coplanar  line  is  transformed 
to  1ft  microstrip  at  the  center  of  the  chip,  and  then  transformed  back  to 
50ft  coplanar  at  the  other  end  of  the  chip.  We  have  thus  far  not  measured 
good  S-parameter  characteristics  with  the  terminated  transformers  and  we 
attribute  this  to  a  problem  in  making  good  contact  with  the  resistors. 

In  fact,  the  measured  dc  resistances  are  extremely  high. 
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One  problem  which  we  discovered  with  the  changeover  from  sapphire 
to  silicon  is  associated  with  holding  the  wafer  during  ion-milling  and 
oxidation.  Since  the  silicon  is  more  fragile  than  sapphire,  one  commonly 
resorts  to  using  vacuum  grease  on  the  back  side  to  hold  the  wafer  to  a 
heat  sink  and  in  t i 1  ted/ i nverted  positions.  We  believe  that  the  cleaning 
steps  in  removing  the  vacuum  grease  caused  contamination  of  the  front 
surface,  particularly  affecting  the  nature  of  contacts  and  junction 
formation.  We  also  suffered  a  failure  in  our  film  thickness  monitor 
which  complicated  our  diagnostic  experiments.  These  problems  appear  to 
have  been  resolved  for  now  in  that  we  are  successfully  fabricating 
junctions  again.  Tests  of  film  contact  resistance,  which  we  believe  is 
the  problem  with  the  terminated  transformer,  have  not  been  completed. 
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6.  PUBLICATIONS 


We  propose  to  submit  a  paper  to  the  1982  Applied  Superconductivity 
Conference  entitled  "SQUID  Vol tage-Control led-Osci 1 lator A  copy  of 
the  proposed  abstract  is  attached  here.  This  paper  will  report  on  the 
analysis  of  the  resistive  SQUID  and  on  the  measurement  which  are 
anticipated  in  the  near  future.  The  Conference  is  November  30  -  December  3, 
1982. 


f  . 
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SQUID  VOLTAGE-CONTROLLED-OSCILLATOR,*  A.H.  Silver, 

R.D.  Sandeil  and  J.Z.  Wilcox,  TRW  —  The  microwave  oscillator 
Is  an  application  of  the  Josephson  junction  which  has  eluded 
successful  development  since  the  earliest  research  on  the 
Josephson  effect.  We  have  investigated  the  SQUID  as  a 
voltage-controlled  source  of  microwave  radiation.  The  low 
impedance  "resistive"  SQUID  can  be  a  relatively  high  power 
(-  nW) ,  tunable  and  monochromatic  source  of  microwaves. 

Studies  of  the  time-dependent  junction  phase  and  the  available 
power  spectrum  as  they  vary  with  such  device  parameters  as  the 
loaded  Q  and  SQUID  8  =  2ttL i  /$<>  not  only  exposes  regions  of 

degenerate  parametric  oscillations  and  chaotic  instabilities, 
but  also  establishes  the  design  rule  for  a  well  behaved 
oscillator,  Q8=2.  Devices  are  fabricated  in  an  integrated 
circuit  format  with  a  monolithic  matching  network  centered  at 
9  GHz.  Methods  of  further  increasing  power  and  impedance  by 
combining  SQUID  have  been  investigated,  particularly  the 
voltage-clamped  dc  SQUID. 


*  Supported  in  part  by  Office  of  Naval  Research  Contract  No. 
N0001^-8l -C-061 5 . 


1.  Category  8. 

2.  Arnold  H.  Silver,  Mail  Stop  RI/1086,  TRW,  One  Space  Park, 
Redondo  Beach,  CA  90278;  213/535-2500. 

3.  A.H.  Silver,  R.D.  Sandeil,  J.Z.  Wilcox. 
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7.  FINANCIAL 


As  of  2  April  1982  we  have  a  total  booked  cost  of  $1*2,273  against 
a  contract  cost  of  $69,055  and  fixed  fee  of  $5,9^5-  We  anticipate 
completing  the  technical  work  by  31  July  with  an  expenditure  of  $67,000. 
$2,055  is  budgeted  for  completing  the  FINAL  REPORT  in  August. 


Figure  1.  Equivalent  circuit  of  the  resistive  SQUID  VCO. 

The  bias  current  is  much  greater  than  the  junction 
current,  1c,  and  the  bias  resistance  is  much 
smaller  ^h»n  the  load  resistance  R. 


Equivalent  circuit  of  the  voltage-clamped  dc  SQUID  with  a  sunning 
output  transformer.  Comparison  with  Fig.  2  shows  that  this  is 
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for  (J”l,  g  =2000  arid  Q=l,5.  This  corresponds  to  w  *2u0.  Time  is 

-1  J 
measured  In  u,  . 
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FIGURE  6.  Computed  power  delivered  to  the  load  resistance  R  in  the 
resistive  SQUIO  at  Uj  as  a  function  of  the  loaded  Q.  Power  is  in  units 
of  (2n)  2L] . 


IGURE  9.  Computed  power  delivered  to  the  load  resistance  In  the  resistive 
SQUfO  with  Q=I  as  a  function  of  g.  Power  is  measured  in  ^L] 

and  frequency  in  to, . 


FIGURE  11.  Composite  of  seven  masks  designed  for  the  four  VCO  circuits. 
Test  junctions  and  resistors  are  placed  above  and  below  each  of  the  2 
VCO's  per  chip. 
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Figure  12.  Photograph  of  transformer  test  chip  mounted  in  a 
brass  holder  and  connected  to  two  OSM  bulkhead  connectors. 
The  tapered  coplanar  lines  at  each  end  are  visible  in  the 
photograph.  The  grey  rectangle  in  the  center  of  the  chip 
is  anodized  Nb  covered  with  SiO.  The  transformers  are  loca¬ 
ted  in  the  same  rectangle  but  are  not  clearly  visible  in  the 
photograph. 


Figure  13.  Response  of  50ft  tapered  coaxial  lines  on  silicon  substrate 
measured  at  4K.  The  upper  graph  shows  S«t  (and  S22 )  f^nt 
curve  is  S,,  at  300K.  The  lower  graph  is  Sot  at  4K. 


